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Abstract 

In the present study, post buckling response of 
hygrothermally induced non-linear elastically supported 
laminated composite plate with uncertain system properties 
using stochastic finite element macromechanical model 
approach is investigated. Uncertain system properties such 
as thermal expansion coefficients, hygroscopic expansion 
coefficients, foundation stiffness parameters, geometric 
properties and lamina material properties are modeled as 
basic random variables. The basic formulation is based on a 
higher-order shear deformation plate theory in von Karman 
nonlinear kinematics. A direct iterative based C°non linear 
finite in conjunction with Taylor series based mean-centered 
first order regular perturbation technique be used to obtain 
mean and standard deviation of post buckling load of 
laminated composite plates. The effects of rising temperature, 
degree of moisture concentrations, boundary conditions, 
plate aspect ratios, total number of plies, nonlinear elastic 
foundation stiffness parameters and fiber orientations are 
analysed. The approach has been validated by comparison of 
results with availability in literature and independent Monte 
Carlo simulation. 
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Introduction 

Nonlinear elastically supported laminated composite 
plates have been used in a great variety of engineering 
applications including aerospace vehicles, aircraft 
runways, launching sites of missile and spacecraft, 
foundation engineering etc. The applications have 
been increased in many folds for the aerospace and 
other engineering fields due to combined features of 



high strength to weight ratio, excellent corrosion 
resistance and good fatigue characteristics. They 
operate in a variety of elevated thermal and moisture 
environments that may have a pronounced impact on 
their performance during their service life. These 
hygrothermal effects are the result of the temperature 
and moisture content variations and are related to the 
difference in the thermal and hygro properties of the 
constituents. The varying environmental conditions 
due to moisture absorption and temperature have an 
adverse effect on the stiffness and strength of the 
composites. As, the matrix is more susceptible to the 
hygrothermal condition than the fiber, the deformation 
is observed to be more in the transverse direction of 
the composite. The rise in moisture and temperature 
reduces the elastic moduli and degrade the strength of 
the materials as well as induces internal initial stresses, 
which may affect the stability as well as the safety of 
the structures. 

As a result, a careful evaluation of environmental 
exposure is required to find the nature and extent of 
their adverse effects upon performance design in 
composite analysis and uncertainties in system 
behavior. Macro-mechanical stochastic modeling 
yields more accurate and exact prediction of system 
behavior and has proved to be superior design. Large 
numbers of parameters are associated with 
manufacturing and fabrication composite components 
as compared to conventional material. For examples 
material and geometric properties, fiber orientations, 
lamina lay-up sequence design and curing parameters. 
Effects of these parameters must be considered in 
accurate prediction of system behavior of composite so 
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that reliability of the structure during it operating life 
can be assessed accurately. Thus there is a need to 
quantify structural system uncertainties in the 
response. This may be appropriately handled by 
modeling system properties in probabilistic sense only. 
Uncertainties in the system properties lead to 
uncertainties in the response behavior of the structure. 
For accurate study of structural behavior, the random 
variation in the system properties should be 
incorporated in the analysis. Otherwise the predicted 
response may differ significantly from the observed 
values and the structures may not be safe. 

In the existing literature of the published works based 
on deterministic analysis using macro mechanical 
model on hygrothermal buckling analysis, notably 
published work of Whitney and Ashton [1971], Flaggs 
et al. [1978], Lee and Yen [1992], Ram and Sinha [1992], 
hygrothermal buckling analysis using micro- 
mechanical model are based on deterministic analysis 
Shen [2001]. Hygrothermal effects on the structural 
behavior of thick laminates applying higher order 
theory of a deterministic analysis has been studied by 
Patel et al. [2002]. In deterministic analysis the system 
parameters are taken to be deterministic and variation 
in the system parameters are ignored which result in a 
conservative design and decrease the potential of the 
composites. It is therefore necessary that the system 
properties are taken as random variables. Limited 
work has been done assuming system properties as 
random and their effects on the performance of 
composites. Nakagiri et al. [1990] studied all edges 
simply supported of a laminate with stochastic finite 
element method taking fiber orientation, layer 
thickness and number of layers as random variables. 
Englested and Reddy [1994] contributed metal matrix 
composites based on probabilistic micro mechanics 
non-linear analysis. Singh at al. [2001] studied the 
buckling of composite cylindrical panels with 
uncertain material properties using probabilistic 
approach. Effects of random system properties on 
initial buckling of composite plates resting on elastic 
foundation using stochastic finite element method was 
studied by Lai et al. [2008]. Post buckling of laminated 
composite plate on elastic foundation with random 
system properties based on HSDT in conjunction with 
FOPT was studied by Singh et al. [2009]. Chen at al. 
[1992] outlined the probabilistic method to evaluate 
the effects of uncertainties in geometric and material 
properties. Relatively little work is available on 
hygrothermal linear buckling analysis of the structures 
made of composites with random system parameters 



using macro-mechanical model Singh and Verma 
[2009]. Pandey et al. [2009] studied the thermo elastic 
stability analysis of laminated composite plates resting 
on nonlinear elastic foundations using analytical 
approach. Singh et al. [2001, 2001] studied the effects of 
random system properties on initial buckling and 
thermal buckling using stochastic finite element 
method. Onkar et al. [2006, 2007] studied the buckling 
problem using stochastic finite element problems. Lai 
et al. [2009] studied the effects on the thermal buckling 
response of laminated composite plates using 
stochastic finite element method, Verma et al.[2009] 
investigated thermal buckling problem with random 
material and geometric properties. 

However, no work dealing with hygrothermal post 
buckling analysis of the laminated composite plate 
with random system properties i.e. material properties, 
geometric properties and nonlinear elastic foundations 
parameters using macro mechanical model is reported 
in the literature to the best of the authors' knowledge. 
The present study is a further extended work of Singh 
et al. [2009] contributions to provide a tool for analysis 
of structural system uncertainties in elevated 
hygrothermal environments especially applicable for 
aerospace engineering and other fields. 

Formulations 

A rectangular arbitrary laminated composite plate 
with of length a, width b, and total thickness h, being 
defined in (X, Y, Z) system with X- and -Y axes located 
in the middle plane and its origin placed at the corner 
of the plate. Let (u, v, w) be the displacement 

paralleling to the (X, Y, Z) respectively as shown in 
Figure 1. The thickness coordinate Z of the top and 
bottom surfaces of any kth layer are denoted by Z(k-i) 
and Z(k) respectively. The fiber of the fah layer is 
oriented with angle 0k to the X- axes. The plate is 
assumed to be attached to the foundation so that no 
separation takes place in the process of deformation. 
The load - displacement relation between the plate 
and the supporting foundations is given as 

P=Kiw+K 3 w 3 -K2V 2 w (1) 

where P is the foundation reaction per unit area, V 2 = 
d 2 / d x 2 + d 2 /dy 2 and Ki, K2 and K3 are linear Winkler 
(normal) foundation, linear Pasternak (shear layer) 
foundation and nonlinear Winkler foundation 
parameters, respectively and w is the transverse 
displacement of the plate. This model is simply known 
as Winkler type when K2 = 0. 
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FIG.l GEOMETRY OF LAMINATED COMPOSITE PLATE WITH 
NONLINEAR ELASTIC FOUNDATION 



Displacement field , Strain Displacement Relations, 
Stress-Strain Relation ,Energy Equations The Strain 
Energy due to Nonlinear Elastic Foundation, Potential 
Energy due to Hygrothermal Stresses, Finite Element 
Formulations are studied in Singh[2009], Reddy[1984], 
Reddy[1997], John[1975], Agrawal[1990], Singh[2002], 
Shankara[1996] and Franklin[1968] . 

The Equation of Motion and its Solution Technique 

Governing Equation, Steps for the Direct Iterative 
Technique and solution Technique are given in Reddy 
[1981], Kleiber [1992], Graham [2001], Yamin [1996], 
Liu [1986] and Graham [2001]. 

Results and Discussions 

The direct iterative based stochastic finite element 
method [DISFEM] is used to illustrate the 
hygrothermally induced post buckling response 
considering different random system parameters of 
laminated composite plate i.e. symmetric and 
antisymmetric. 

In this section the mean post buckling load and effect 
of randomness in the system properties on the 
postbuckling load of geometrically non linear 
laminated composite plate resting on nonlinear elastic 
foundations are investigated in hygrothermal 
environments. The present methodology for mean and 
coefficient of variation of post buckling load in 
hygrothermal environment is verified with available 
results in literature and independent MCS. The lamina 
coefficients of hygroscopic expansion, thermal 
expansion and foundation stiffness parameters 
including geometric and material properties are 
modeled as basic random variables (RVs) as stated 
earlier. However, the results are only presented taking 
coefficient of variation (COV) of the system property 
equal to 0.10 Liu [1986] as the nature of the SD 
(standard deviation) variation is linear and passes 
through the origin. Hence, the presented results would 



be sufficient to extrapolate the results for other COV 
value keeping in mind the limitation of the FOPT 
Singh et al. [2009,2002]. The thickness of all the lamina 
is assumed to be constant and of same material 
without varying individual property of materials used. 
The obtained results have been compared with MCS 
and those available in the literature. For the present 
study a nine nodded serendipity element with 63 
degree of freedom (DOFs) per element, a modified 
form HSDT based C° finite element model has been 
used for discretizing the laminate. Based on 
convergence, an (8x8) mesh has been used throughout 
for evaluation of the results. Programme in Mat lab is 
developed for the computational analysis. Results 
have been computed by employing the full (3x3) 
integration rule for bending stiffness matrices, 
hygrothermal load vector and the geometric stiffness 
matrices and the reduced (2x2) integration rule for 
computing the shear stiffness matrices to avoid shear 
locking in the thin plates. 

All edges simply supported (SSSS (SI)): 

u = w = 6 = y/ = 0, at x = 0, a; v = w = = y/ x = at y = 0, b ; 

All edges simply supported (SSSS (S2)): 

v = w = 9 = iff = 0, at x = 0, a; u = w = 6 X =y/ t = at y = 0, b ; 

All edges clamped (CCCC): 

u = v = w = y/ x =y/ =0 y = 9 , = 0, at x = 0, a and y = 0, b; 

The used plate geometry is characterized by aspect 
ratios (a/b) = land 2, side to thickness ratios (a/h) = 10, 
20, 30, 40, 50, 60, 80 and 100. The mean values of 
coefficients of hygroscopic expansion, coefficients of 
thermal expansion, geometric properties, material 
constants and foundation stiffness parameters are used 
for computation. We consider now a second steps as 
the elastic constants, thermal expansion coefficients 
and coefficients of hygroscopic expansion of each 
layers are assumed to be linear function of 
temperature and moisture. The only exception is the 
Poisson's ratio, which can be assumed as constant due 
to weakly dependency on temperature and moisture 
change. The coefficient of hygroscopic expansion jSnin 
longitudinal direction is considered to be zero because 
moisture contents have negligible effects on fiber and 
have effects on matrix in transverse direction. For the 
coefficient of variation of material and geometric 
properties varying from to 20% with a/h= 10, 20, 30, 
40, 50, 80 and 100 subjected to uniform temperature 
and moisture rise have been considered. The lamina 
material properties such as En, E22, G12, G13, G23, m, az, 
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02 and h are modeled as basic RVs. En, and £22 are 
longitudinal and transverse elastic modulii, 
respectively, G12, G13, G23 are shear modulii, m and «2 
are longitudinal and transverse coefficients of thermal 
expansion, f>\ and 02 are longitudinal and transverse 
coefficients of hygroscopic expansion and h are the 
total thickness of the composite laminate considered 
for the analysis. The input random variables fa is 
related and sequenced as 

b\= En, t>2=£22, fo=Gi2, bi= G13, fc=G23, be=ai, bi=ai, bs= p2 
and£>9=/z. The input random variables for plate resting 
on nonlinear elastic foundations are related as bi= En, 
bi=E2-i, bs=Gn, bi= G13, bs=Gi3, bs=a\, bi=ai, bs= 02, b<s=k\, 
bw=ki and bn=ki. The non dimensionalised foundation 
parameters used for the analysis are as 

k, = K 3 a 4 1 E' 2 ' 2 h 3 ; k 2 = K 2 a 4 1 E^ 2 h 3 ; 

k 3 = K 3 a 4 I E d 27 h 2 . 

Table 1(a), 1(b) and 1 (c) are compared for validation of 
mean post buckling analysis. Tables 2(a), 3(a), 4(a), 5(a) 
and 7(a) give generated results for various 
combinations of random variables with geometric 
parameters for second order statistics for post buckling 
load of laminated composite plate in hygrothermal 
environment. The numerical values and relationship 
between the mean values of the material properties for 
graphite epoxy composite with change in temperature 
and moisture. 

For plates resting on nonlinear elastic foundations 
Tables 1(d), 1(e) and 1(f) are compared for validation of 
mean post buckling analysis. Tables 2(b), 3(b), 4(b), 5(b), 
5(c), 6(a), 6(i), 7(b) and 7(c) give generated results for 
various combinations of random variables with 
geometric parameters for second order statistics for 
post buckling load of laminated composite plate 
resting on nonlinear elastic foundations in 
hygrothermal environment. The numerical values and 
relationship between the mean values of the material 
properties for graphite epoxy composite with change 
in temperature and moisture. 

Validation Study: Mean and Standard Deviations 
Without Foundations 

In order to validate the proposed outlined approach, 
the results for the mean and standard deviation are 
compared with those available in the literature and an 
independent Monte Carlo simulation technique. 

1) Mean Hygrothermal Post Buckling Analysis 

Table 1 (a) & Table 1(b) show the comparison of 
effect of temperature and moisture respectively on 



the critical load Acri= Nxcr/ (Nxcr)c=o% or t=3ook of anti- 
symmetric cross ply (0°/90°)2t, clamp supported 
laminate, side to thickness ratio a/h =10 with those 
available results. The results are in good agreement. 
Table 1(c) shows the comparison of effect of 
moisture and temperature on the dimensionalised 
buckling load of cross ply [0°/90°], clamped 
supported square laminate with those available 
results of Singh and Verma.[2009]. The present 
results are in close tolerance with the available 
literature. Table 1 (d), 1(e) show the comparison of 
effect of temperature and moisture respectively on 
the critical load Acri= Nxcr/ (Nxcr)c=o% or t=3ook of anti- 
symmetric cross ply (0°/90°)2t, simply supported S2 
laminate, plate thickness ratio (a/h) =10 with those 
available results SaiRam and Sinha [1992]. The 
results are in good agreement. Table 1(f) shows the 
comparison of parametric studies results for aspect 
ratio a/b=l, thickness ratios a/h=10, simple support 
SSSS (S2), biaxial compression , angle ply [±45°]2t, 
on dimensionless mean thermal post buckling load 
of laminated composite plate subjected to 
temperature and moisture independent (TID) 
material properties. It is observed that foundation 
parameters have significant effects as shown by 
bounds of buckling load and post buckling strength. 
The material properties of the composite plate are: 
E1/E2 =25; Gi2=0.5E2; G23=0.2E 2 ; E2o=l*le5; zm=0.25; 
Eio=mi*E2o; Gi2o=m2*E2o; Gi3o=Gi2o; G23o=m3*E2o; 
P2i=in2*E2o/£io; ano =l*le-6; ana =10.0*le-6; an = ; 
ao =le-6; The non dimensionalised foundation 
parameters are as; Ki= fci*Dii / a 4 , Ki= fo*Dn / a 2 , K3= 
fo*Dn / a 4 *h 2 and (Tcr = Ti* ao*1000). The results 
obtained are in good agreements and validated 
with the available literature results of Pandey et al. 
[2009]. 

2) Validation Result for Random Material 
Properties 

As mentioned earlier, no results are available in 
reported literature for variance of post buckling 
load in hygrothermal environment using HSDT 
based DISFEM. The verification of the results for 
COV of hygrothermal post buckling load is 
compared with available results obtained by MCS 
approach using samples of random number 
generation. The present direct iterative method in 
conjunction with first order perturbation technique 
[DISFEM] results for random response has been 
validated with an independent MCS approach. 
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alb 




Temperature T(K) 


300 


325 


350 


375 


400 


425 


0.5 


SaiRam and Sinha 
[1992] 


1.000 


0.988 


0.951 


0.911 


0.888 


0.865 


Present 


1.000 


0.991 


0.954 


0.913 


0.890 


0.868 


2.0 


SaiRam and Sinha 
[1992] 


1.000 


0.994 


0.931 


0.865 


0.828 


0.792 


Present 


1.000 


0.995 


0.930 


0.862 


0.827 


0.792 


a/b 




Moisture Concentration C(%) 


0.00 


0.25 


0.50 


0.75 


1.00 


1.25 


1.50 


0.5 


^aiT?3m an/i Qinna 

iJJClllVClll I Cli Lt_l iJJll 11 let 

[1992] 


1.000 


0.992 


0.984 


0.977 


0.970 


0.965 


0.960 


Present 


1.000 


0.9916 


0.9833 


0.9752 


0.9672 


0.960 


0.953 


2.0 


SaiRam and Sinha 
[1992] 


1.000 


0.997 


0.993 


0.989 


0.987 


0.985 


0.982 


Present 


1.000 


0.997 


0.994 


0.991 


0.989 


0.987 


0.985 



Table 1(c) 



a/h 


Singh and Verma [2009] 


Present [HSDT] 


Moisture concentration C 
(%) 


Temperature T 
(K) 


buckling load Nx 
(N/m) (10 8 ) 


Ei & G12 


buckling load Nx (N/m) x 
(10 8 ) 


10 


0.1 


300 


7.475 


9.5 /6.0 


7.4412 


20 


0.5 


350 


1.1536 


9.0/6.0 


1.5837 


50 


1.0 


400 


0.097 


8.5/ 6.0 


0.1102 



Table 1 (d) 



a/b 




Temperature T(K) 


300 


325 


350 


375 


400 


425 


0.5 


SaiRam and 
Sinha [1992] 


1.000 


0.983 


0.955 


0.926 


0.907 


0.889 


Present 


1.000 


0.986 


0.963 


0.939 


0.923 


0.910 


Table 1 (e) 



a/b 




Moisture Concentration C(%) 


0.00 


0.25 


0.50 


0.75 


1.00 


1.25 


1.50 


0.5 


SaiRam and 
Sinha [1992] 


1.000 


0.982 


0.965 


0.949 


0.934 


0.920 


0.906 


Present 


1.000 


0.984 


0.968 


0.953 


0.939 


0.926 


0.912 



Table 1(f) 



(Wmax/fe) 


Pandey et al.[2009] 


Present[HSDT] 


0.0 


9.850 


9.981 


0.2 


10.076 


10.157 


0.4 


10.998 


11.550 


0.6 


12.282 


11.898 
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0.06 
0.07 
0.0E 

_0.0£ 

J 



a o.os 



0.01 
0. 



— ■— : Present [DISFEM] 
Present [MCS] 

(a)h/h=10, a/b=1, 
SSSS(S2),Wmax/h=0.4 

[0°/90°] 2T , TIDcase 

^,=0.6684x1 8 N/m 

\ Jn |=0.7643x1 8 N/m 
T=325K, C=0.010 




"OTTO 0TT5 - 

SD/mean, 



0.08 
0.07 
0.06 
_ 0.05 
... 0.04 
0.03 
0.02 
0.01 



— ■— : Present [DISFEM] 
— •— : Present [MCS] 

(b) Ml=10, a/b-1, SSSS(S2), (T.T) 
[0°/90 /0%°] 
X crr 4.1542x10 8 W/m 
/. crn p4.2514x10 8 N/m 
T-325K, C-0.25% 



0.10 0.15 
SD/mean, h 



Figure. 2 Validation of present DISFEM results with 
independent MCS results for: (a) Random material 
properties (b) Random geometric properties. 

Fig. 2 (a) plots the normalized standard deviation, 
coefficient of variation of hygrothermal post 
buckling load versus the coefficient of variation of 
the random material constants of cross-ply 
[0 /90°/0°/90 ] square laminated composite plates, 
SSSS (S2), a/h=20 with amplitude ratio (Wmax/h =0.2) 
subjected to combination of uniform temperature 
and transverse temperature distribution changing 
from to 20%. It is assumed that one of the material 
properties (i.e., En) changes at a time keeping other 
as a deterministic, with their mean values of the 
material properties. The dashed line is the present 
[DISFEM] result and the solid line is independent 
MCS approach. For the MCS approach, the samples 
are generated using Mat Lab to fit the mean and 
standard deviation. These samples are used in 
response equation Eq. (25), which is solved 
repeatedly, adopting conventional eigenvalue 
procedure, to generate a sample of the 
hygrothermal post bucking load. The number of 
samples used for MCS approach is 10,000 based on 



satisfactory convergence of the results. The normal 
distribution has been assumed for random number 
generations in MCS. It can also be observed that the 
present DISFEM results are in good agreement with 
MCS result. 

3) Validation for Random Geometric Properties 

Figure 2(b) demonstrates a comparison of the 
obtained result in the present DISFEM and 
independent MCS approach treating random 
lamina thickness as random variables keeping 
others as constant. The composite laminate having 
stacking sequence [0°/90°]2t thickness ratio a/h= 10 
and amplitude ratio Wmax/h=0.2 are used in the 
analysis. For the range of COV in the geometric 
properties such as lamina plate thickness (h) 
changing at a time keeping other as deterministic 
i.e. mean values are considered i.e. to 20%. It is 
concluded that DISFEM as shown by dotted line 
gives acceptable results with the solid line of 
independent MCS using 12,000 samples of random 
number generation for the range of dispersion 
taken in the study. 

Parametric Analysis of Second Order Statistics 
Without Foundations 

Table 2 shows the effects of amplitude ratios (Wmax/h) 
and nonlinear elastic foundation parameters with 
individual random system property h, {(i =1 to 11) = 
0.10} on the dimensionlised mean (T*cmi = Tcmi e+08) 
N/m and dimensionless coefficient of variation (A.cmi) of 
hygrothermal post buckling load of cross ply [0°/90°]2t 
square laminated composite plates, subjected to 
uniform constant temperature and moisture (U.T) and 
in-plane uni-axial compression, plate thickness ratio 
(a/h) =20, rise in temperature and moisture AT= 325K, 
AC= 0.010, with simple support S2 boundary 
conditions. (T*cri = Tcrie+08) - mean dimensionalized 
linear solution. It is seen that the effects of amplitude 
ratios increases the mean values and COV for plate 
resting on nonlinear Winkler foundation. These values 
further increase when the plate is resting on Pasternak 
nonlinear elastic foundation. The mean values and 
COV are higher for input random variable En for plate 
resting on (fa=100, fa=10, fe=200) Pasternak nonlinear 
elastic foundation. However the effects on coefficient 
of thermal expansion an are fever compared to all 
other input random variables . 

Table 3(a) shows the effects of plate thickness ratios 
(a/h) , amplitude ratios (Wmax/h) and the random input 
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variables h, {i =(1..8), (6, 7), (8) and (9) = 0.10}.on the 
dimensionlised mean (T*cmi = Tcmie+08) N/m and 
dimensionless coefficient of variation (kani) of 
hygrothermal post buckling load of cross ply [0°/90°]2t 
square laminated composite plates with simple 
support (SI) and clamped support boundary 
conditions, rise in temperature and moisture (AT=325K, 
AC=0.25%). It is seen that for simple support SI 
conditions, the mean hygrothermal post buckling load 
decreases as the side to thickness ratio increases while 
on increasing the amplitude ratio the result is reversed. 
The effects of combined random system property 
except lamina thickness with increase side to thickness 
ratio and amplitude ratio the coefficient of variation of 
hygrothermal post buckling load of the plate increases. 
In the case of plate with clamped support conditions, 
the mean value with different combinations increase 
much more compared to plate with simple support 
conditions. However, the effects on COV with 
different random system properties vary significantly 
in case of plate with clamped support conditions 
keeping same temperature and moisture conditions. 

Table 3(b) shows the effects of plate thickness ratios 
(a/h), nonlinear elastic foundation parameters and 
amplitude ratios (Wmax/h) with random input variables 
h,[{(i =1.. . 11), (6, 7),(8) and (9..11)} = 0.10] on the 
dimensionalised mean (T*cmi = Tcmi e+08) N/m and 
dimensionless coefficient of variation (lemi) of 
hygrothermal post buckling load of cross ply [0°/90°]2t , 
square laminated composite plate, subjected to 
uniform constant temperature and moisture , in-plane 
uni-axial compression with simple support S2 
boundary conditions. AT= 325K, AC= 0.010, (T*cri = 
Tcrie+08) N/m - mean dimensionalised linear solution 
in brackets. It is observed that for a given plate 
thickness ratio the effects of increase in amplitude ratio 
increases the mean values for the plate resting on 
nonlinear elastic foundations and these values further 
increase for Pasternak nonlinear elastic foundations. 
The effects on COV are of significance for all combined 
random input variables. If the plate thickness ratio is 
further increased then mean values further decreases 
and COV is significant for all combined random input 
variables and Winkler nonlinear elastic foundations 
compared to Pasternak nonlinear elastic foundations. 

Table 4 shows the effects of aspect ratios (a/b) , 
nonlinear foundation parameters and amplitude ratios 
(Wmax//i) with random input variables bi,[{(i =1 .. 11), (6, 
7),(8) and (9..11)} = 0.10 on the mean (T*cmi = Tcmi e+08) 
N/m and dimensionless coefficient of variation (tami) of 
hygrothermal post buckling load of angle ply [±45°]2t 



laminated composite plate, subjected to uniform 
constant temperature and moisture, in-plane uni-axial 
compression with a/h=30. AT= 325K, AC= 0.010. (T*cri = 
Tcrie+08) - mean linear solution. It is noticed that for 
square plate the effects of amplitude ratios increases 
the mean values for plate resting on Winkler nonlinear 
elastic foundations and these values are further higher 
for plate resting on Pasternak nonlinear elastic 
foundations, however there is appreciable increase in 
coefficient of variations. For rectangular plate the 
mean values sharply decline and COV increase 
significantly and is dominant for all combined input 
random variables. The COV is also higher for f>27. when 
plate rests on Winkler nonlinear elastic foundations . 

Table 5(a). shows the effects of support conditions, 
nonlinear foundation parameters and amplitude ratios 
(Wmax//i) with random input variables bi,[{(i =1 .. 11), (6, 
7),(8) and (9. .11)} = 0.10] on the dimensionlised mean 
(T*cmi = Tcmi e+08) N/m and dimensionless coefficient of 
variation (Xani) of hygrothermal post buckling 
load, ,plate thickness ratio (a/h) =20, of angle ply 
[±45°]2T laminated square composite plate, subjected to 
of uniform constant temperature and moisture, in- 
plane uni-axial compression, AT= 325K, AC= 0.010, 
(T*cri = Tcrie+08) - mean linear solution. It is seen that for 
simple support SSSS (SI) square composite plate 
resting on Winkler nonlinear elastic foundation the 
increase in amplitude ratio increases the mean values 
and COV decreases. Mean values are further higher for 
plate resting on Pasternak nonlinear elastic 
foundations. The combined effects of all input random 
variables are dominant for Pasternak nonlinear elastic 
foundations. When the plate is simply supported SSSS 
(S2) and resting on Winkler or Pasternak nonlinear 
elastic foundations there is marginal variation in mean 
values and COV. 

Table 5(b) shows the effects of support conditions, 
nonlinear foundation parameters and amplitude ratios 
(Wmax//j) with random input variables bi,[{(i =1 .. 11), (6, 
7),(8) and (9..11)} = 0.10] on the mean (T*cmi = Tcmi e+08) 
N/m and dimensionless coefficient of variation (kcmi) of 
hygrothermal post buckling load, ,plate thickness ratio 
(a/h) =20, of angle ply [±45°]2t laminated square 
composite plate, subjected to of uniform constant 
temperature and moisture, in-plane uni-axial 
compression, AT= 325K, AC= 0.010, (T*cri = Tcrie+08) - 
mean linear solution. It is seen that for clamped 
support CCCC square composite plate resting on 
Winkler nonlinear elastic foundation the increase in 
amplitude ratio significantly increase the mean values 
and COV decreases. Mean values are further higher for 
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TABLE 2 







COV, Ui 




Wmax/fc 


(fcWOO, ki=00 ,fo=100) 


(fa=100, fc=10, fe=100) 


(fci=100, fc=10, fo=200) 


En (i=l) 


0.2 


(Tcmi = 0.1887) 0.1053 


(Tcmi =0.3324) 0.1523 


(Tcmi =0.3345) 0.1550 




0.4 


(Tc,„i = 0.2057) 0.1023 


(Tcmi =0.3803) 0.1482 


(Tcmi =0.3881)0.7227 




0.6 


(Tcmi = 0.2290) 0.0965 


(Tcmi =0.4128)0.2339 


(Tcmi =0.4947)0.2274 




Tcrl 


(0.1826) 


(0.3003) 


(0.3003) 


E22 (i=2) 


0.2 


0.0677 


0.0589 


0.0586 




0.4 


0.0644 


0.0594 


0.0581 




0.6 


0.0633 


0.0533 


0.0579 


Gi2(/=3) 


0.2 


0.0098 


0.0052 


0.0052 




0.4 


0.0091 


0.0050 


0.0049 




0.6 


0.0081 


0.0062 


0.0057 


G13 (i=4) 


0.2 


0.0023 


0.0068 


0.0068 




0.4 


0.0024 


0.0077 


0.0076 




0.6 


0.0036 


0.0048 


0.0100 


G23 (!=5) 


0.2 


0.0012 


0.0034 


0.0034 




0.4 


0.0012 


0.0038 


0.0038 




0.6 


0.0018 


0.0024 


0.0050 




0.2 


0.0014 


7.98e-04 


7.93e-04 


an (1=6) 


0.4 


0.0013 


6.97e-04 


6.83e-04 




0.6 


0.0012 


6.42e-04 


5.36e-04 


a22 (i'=7) 


0.2 


0.0110 


0.0062 


0.0062 




0.4 


0.0101 


0.0054 


0.0053 




0.6 


0.0090 


0.0050 


0.0042 


/?22 (i=8) 


0.2 


0.0635 


0.0360 


0.0358 




0.4 


0.0582 


0.0315 


0.0309 




0.6 


0.0523 


0.0290 


0.0242 


fci (»=9) 


0.2 


0.0316 


0.0047 


0.0045 




0.4 


0.0283 


0.0039 


0.0038 




0.6 


0.0221 


0.0120 


0.0027 


fa (WO) 


0.2 





0.0225 


0.0222 




0.4 





0.0196 


0.0191 




0.6 





0.0263 


0.0150 


fo(Wl) 


0.2 


7.13e-04 


0.0013 


0.0024 




0.4 


0.0026 


0.0042 


0.0078 




0.6 


0.0044 


0.0167 


0.0338 



plate resting on Pasternak nonlinear elastic 
foundations. The combined effects of all input random 
variables are dominant for Pasternak nonlinear elastic 
foundations. When the plate is CSCS supported and 
resting on Winkler or Pasternak nonlinear elastic 
foundations there is decrease in mean values and 
increase in COV. The influence of Clamped support 
with the plate resting on Pasternak nonlinear elastic 
foundations is dominant. 



Table 6(a) shows the effects of no of layers, nonlinear 
foundation parameters and amplitude ratios (Wmax/ft) 
with random input variables bi,[{(i =1 .. 11), (6, 7),(8) 
and (9. .11)} = 0.10] on the dimensionlised mean (T* C mi = 
Tcmi e+08) N/m and dimensionless coefficient of 
variation (Xani) of hygrothermal post buckling load, 
plate thickness ratio (a/h) =40, of angle ply laminated 
square composite plate subjected to simple boundary 
support S2, uniform constant temperature and 
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moisture, in-plane uni-axial compression, rise in 
temperature and moisture AT= 325K, AC= 0.010, (T*cri = 

Tcrie+08) - mean dimensionalized linear solution. It is 
observed that for a lay-up cross ply laminated 
composite plate simply supported and resting on 
Winkler nonlinear elastic foundations in study, the 
increase in amplitude ratio increases the mean values 
and COV decreases. The mean values are further 
increased and COV is decreased when the plate is 
resting on Pasternak nonlinear elastic foundations. The 
combined effects of all random variables are dominant 
besides ^22 effects on coefficient of variations. With the 
increase of number of layers, mean values further 
increases marginally and COV decrease. The mean 
values and COV are higher compared to cross ply 
laminated composite plate with similar foundation 
parameters. 

Table 6 (b) shows the effects of no of layers, nonlinear 
foundation parameters and amplitude ratios (Wmax//i) 
with random input variables bi,[{(i =1 •• 11)/ (6, 7),(8) 



and (9..11)} = 0.10] on the mean (T*cmi = Tcmi e+08) N/m 
and dimensionless coefficient of variation (Xcmi) of 
hygrothermal post buckling load, plate thickness ratio 
(a/h) =40, of angle ply laminated square composite 
plate subjected to simple boundary support S2, 
uniform constant temperature and moisture, in-plane 
uni-axial compression, rise in temperature and 
moisture AT= 325K, AC= 0.010, (T*cri = Tcrie+08) - mean 
linear solution. It is noticed that for angle ply 
laminated composite plate simply supported and 
resting on Winkler nonlinear elastic foundations in 
study, the increase in amplitude ratio increase the 
mean values and COV decrease. The expected mean 
values are further increased and coefficient of 
variations is decreased when the plate is resting on 
Pasternak nonlinear elastic foundations. The combined 
effects of all random input variables are dominant 
besides ^22 effects on coefficient of variations. With the 
increase of number of layers mean values further 
increase marginally. 



TABLE 3(A) 











COV, Xcml 


BCs 


a/h 


Wmax/h 


(Tend) 


bi 








(i=1..8) 


(i=6,7) 


(i=8) 


(i=9) 






0.2 


0.0112 


0.2040 


0.0664 


0.0922 


0.0353 




40 


0.4 


0.0124 


0.1895 


0.0600 


0.0833 


0.0326 






0.6 


0.0141 


0.1739 


0.0526 


0.0730 


0.0295 


SSSS(Sl) 




0.2 


0.0034 


0.3014 


0.1474 


0.2047 


0.0354 




60 


0.4 


0.0037 


0.2755 


0.1331 


0.1848 


0.0324 






0.6 


0.0042 


0.2462 


0.1165 


0.1618 


0.0290 






0.2 


0.0014 


0.4694 


0.2576 


0.3578 


0.0349 




80 


0.4 


0.0016 


0.4271 


0.2332 


0.3239 


0.0320 






0.6 


0.0018 


0.3777 


0.2046 


0.2841 


0.0286 






0.2 


0.0356 


0.2481 


0.0208 


0.0290 


0.0403 




40 


0.4 


0.0376 


0.2407 


0.0197 


0.0274 


0.0392 






0.6 


0.0407 


0.2319 


0.0182 


0.0253 


0.0379 


cccc 




0.2 


0.0113 


0.2405 


0.0439 


0.0610 


0.0412 




60 


0.4 


0.0119 


0.2323 


0.0415 


0.0577 


0.0396 






0.6 


0.0130 


0.2218 


0.0382 


0.0530 


0.0374 






0.2 


0.0049 


0.2524 


0.0758 


0.1052 


0.0419 




80 


0.4 


0.0052 


0.2424 


0.0715 


0.0994 


0.0400 






0.6 


0.0056 


0.2290 


0.0656 


0.0911 


0.0373 
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TABLE 3(B) 











a/h=l0 










a/h=30 






Foundation 


Wmax/fo 






















parameters 




















































COV, A.cml 










COV, Xcm\ 








(Tcrnl) 


bi 


(Tend) 


bi 








(i=l..ll) 


(i-6,7) 


0=8) 


(i=9,ll) 




(i=l..ll) 


(i-6,7) 


(i=8) 


(i=9,ll) 


(h=100, 


0.2 


0.7673 


r\ *i t r\r\ 

O.llOO 


0.0027 


0.0156 


0.0087 


0.1170 


0.1826 


0.0179 


0.1024 


0.0238 


/C2— UU ,/C3— 


0.4 


0.8318 


0.2065 


0.0025 


0.0144 


0.0304 


0.1464 


0.1514 


0.0143 


0.0818 


0.0088 




0.6 


0.8689 


0.1609 


0.0024 


0.0138 


0.0211 


0.1690 


0.1508 


0.0124 


0.0709 


0.0110 




Tcrl 


(0.6684) 










(0.0835) 










(fa=100, fc=10. 


0.2 


1.0354 


0.H53 


0.0020 


0.0231 


0.0299 


0.1583 


0.1899 


0.0132 


0.0757 


0.0233 


/C3— 


0.4 


1.2473 


0.H57 


0.0021 


0.0119 


0.0299 


0.1819 


0.1763 


0.0115 


0.0658 


0.0196 




0.6 


1.1632 


0.1088 


0.0017 


0.0103 


0.0284 


0.2069 


0.1308 


0.0101 


0.0579 


0.0183 




Tcrl 


(0.9855) 










(0.1412) 










(fa=100, fo=10. 


0.2 


1.0664 


0.H28 


0.0020 


0.0112 


0.0288 


0.1597 


0.1436 


0.0131 


0.0750 


0.0227 


fc>=200) 


0.4 


1.1353 


0.1166 


0.0018 


0.0105 


0.0308 


0.1859 


0.1370 


0.0112 


0.0644 


0.0202 




0.6 


1.0489 


0.H99 


0.0020 


0.0114 


0.0319 


0.2134 


0.1374 


0.0098 


0.0561 


0.0216 




Tcrl 


(0.9855) 










(0.1359) 











TABLE 4 







a/b=l 


a/b=2 










COV, Xcml 








COV, Xcml 




Foundation 
parameters 


Wmax/ft 






















(Tcrnl) 


bi 


(Tcrnl) 


bi 








(HL..11) 


0=6,7) 


0=8) 


(i=9..11) 




(i=l..ll) 


0=6,7) 


0=8) 


(i=9..11) 


(fa=100, 


0.2 


0.1329 


0.1503 


0.0157 


0.0901 


0.0051 


0.0463 


0.2804 


0.0451 


0.2587 


0.0022 


fa=00 fo=100) 


0.4 


0.1555 


0.1385 


0.0134 


0.0770 


0.0040 


0.0684 


0.2012 


0.0305 


0.1749 


0.0019 




0.6 


0.1775 


0.1355 


0.0118 


0.0674 


0.0055 


0.0835 


0.1766 


0.0250 


0.1434 


0.0019 




Tcrl 


(0.1107) 










(0.0339) 










(fa=100, 


0.2 


0.1521 


0.1584 


0.0137 


0.0787 


0.0138 


0.0512 


0.2637 


0.0408 


0.2337 


0.0099 


fa=10, 
fa=100) 


0.4 


0.1680 


0.1499 


0.0124 


0.0713 


0.0128 


0.0732 


0.1000 


0.0285 


0.1635 


0.0067 




0.6 


0.1891 


0.1411 


0.0110 


0.0633 


0.0115 


0.0881 


0.1749 


0.0237 


0.1359 


0.0055 




Tcrl 


( 0.1413) 










(0.0389) 










(fa=100. 


0.2 


0.1524 


0.1589 


0.0137 


0.0785 


0.0138 


0.0514 


0.2636 


0.0406 


0.2329 


0.0099 


h=W, 
fa=200) 


0.4 


0.1692 


0.1520 


0.0123 


0.0708 


0.0129 


0.0737 


0.1972 


0.0283 


0.1624 


0.0071 




0.6 


0.1917 


0.1453 


0.0109 


0.0625 


0.0122 


0.0888 


0.1760 


0.0235 


0.1348 


0.0061 




Tcrl 


(0.1413) 










(0.0389) 
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TABLE 5 (A) 







SSSS (SI) 


SSSS (S2) 






(Tcrnl) 




COV, Xcml 




(Tcrnl) 




COV, A,c,„i 




Foundation 


WmWIl 










































parameters 








bi 








bi 










(i=l..ll) 






(i=9..11) 




(i=l..ll) 


<i=6 71 




(i=9..11) 


(fa=100. 


0.2 


0.2794 


0.1280 


0.0075 


0.0429 


0.0044 


0.2744 


0.1302 


0.0076 


0.0436 


0.0037 


fc=00 ,fe=100) 


0.4 


0.3221 


0.1236 


0.0065 


0.0372 


0.0042 


0.3176 


0.1255 


0.0066 


0.0377 


0.0040 




0.6 


0.3438 


0.1359 


0.0061 


0.0348 


0.0101 


0.4141 


0.1507 


0.0050 


0.0289 


0.0139 




Tcrl 


(0.2387) 










(0.2381) 










(fa=100, 


0.2 


0.3369 


0.1580 


0.0062 


0.0355 


0.0199 


0.4890 


0.1313 


0.0043 


0.0245 


0.0126 


fa=10, fa=100) 


0.4 


0.3471 


0.1402 


0.0060 


0.0345 


0.0146 


0.3915 


0.1464 


0.0053 


0.0306 


0.0163 




0.6 


0.4023 


0.1373 


0.0052 


0.0298 


0.0141 


0.3376 


0.1476 


0.0062 


0.0355 


0.0169 




Tcrl 


(0.3033) 










(0.2982) 










(fa=100, 


0.2 


0.3490 


0.1384 


0.0060 


0.0343 


0.0147 


0.4368 


0.1245 


0.0048 


0.0274 


0.1115 


fo=10, fa=200) 


0.4 


0.4221 


0.1328 


0.0049 


0.0284 


0.0124 


0.3267 


0.1543 


0.0064 


0.0366 


0.0183 




0.6 


0.4081 


0.1397 


0.0051 


.0293 


0.0119 


0.3281 


0.1520 


0.0064 


0.0365 


0.0158 




Tcrl 


(0.3033) 










(0.2982) 











TABLE 5 (B) 







cccc 


CSCS 






(Tcrnl) 




COV, Xcml 




(Tcrnl) 




COV, A,cml 




Foundation 












































parameters 








bi 








bi 










(i=l..ll) 


(i=6,7) 


(i=8) 


(i-9..11) 




(i=l..ll) 


(i-6,7) 


(i=8) 


(i=9..11) 


(fci=100. 


0.2 


0.4227 


0.1285 


0.0049 


0.0283 


0.0017 


0.3632 


0.1331 


0.0058 


0.0330 


0.0051 


fa=00 ,fo=100) 


0.4 


0.5262 


0.1331 


0.0040 


0.0228 


0.0038 


0.3256 


0.1359 


0.0064 


0.0368 


0.0049 




0.6 


0.5130 


0.1356 


0.0041 


0.0233 


0.0055 


0.4058 


0.1325 


0.0051 


0.0295 


0.0054 




Tcrl 


(0.3758) 










(0.2878) 










(fci=100. 


0.2 


0.4691 


0.1354 


0.0045 


0.0511 


0.0096 


0.4114 


0.4114 


0.0051 


0.0291 


0.0142 


fo=10, fa=100) 


0.4 


0.5375 


0.1324 


0.0039 


0.0223 


0.0089 


0.3753 


0.1480 


0.0056 


0.1319 


0.0140 




0.6 


0.4938 


0.1361 


0.0042 


0.0242 


0.0098 


0.4591 


0.1427 


0.0045 


0.0261 


0.0127 




Tcrl 


(0.4306) 










(0.3433) 










(fci=100. 


0.2 


0.5521 


0.1359 


0.0038 


0.0179 


0.0095 


0.4124 


0.1467 


0.0051 


0.0290 


0.0140 


fc=10, fe=200) 


0.4 


0.5466 


0.1300 


0.0038 


0.0219 


0.0086 


0.3777 


0.1489 


0.0055 


0.0317 


0.0140 




0.6 


0.5338 


0.1313 


0.0039 


0.0224 


0.0096 


0.4064 


0.1356 


0.0051 


0.0295 


0.0115 




Tcrl 


(0.4306) 










(0.3433) 
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Table 7(a) shows the effects of environmental 
conditions, nonlinear foundation parameters and 
amplitude ratios (Wmax/h) with random input variables 
h,[{(i =1 •• 11), (6, 7),(8) and (9..11)} = 0.10] on the 
dimensionlised mean (T*cmi =Tcmie+08) N/m and 
dimensionless coefficient of variation (Xcmi) of 
hygrothermal post buckling load of angle ply [±45°]2t 
laminated square composite plates with SSSS (S2) 
boundary conditions subjected to in-plane uni-axial 
compression , plate thickness ratio (a/h)=20. (T*ai = 
Tcrie+08) - mean dimensionalized linear solution. It is 
seen that for given environmental conditions the 
increase in amplitude ratios increases the mean values 
and coefficient of variations decreases. The combined 
effects of all random input variables increase the 
coefficient of variations when plate is resting on 
Winkler nonlinear elastic foundations. On change of 
foundation parameters to Pasternak these values 
further increase and coefficient of variations decreases. 
With the increase of moisture contents keeping 
temperature constant the effects of amplitude ratios 
decrease the mean values for Winkler nonlinear elastic 
foundations, however these values are higher for 
Pasternak nonlinear elastic foundations. 

Table 7(b) shows the effects of environmental 



conditions, nonlinear foundation parameters and 
amplitude ratios (Wmax/h) with random input 
variables h,[{(i =1 .. 11), (6, 7),(8) and (9..11)} = 0.10] 
on the dimensionlised mean (T*cmi =Tcmie+08) N/m 
and dimensionless coefficient of variation (Xcmi) of 
hygrothermal post buckling load of angle ply [±45°]2t 
laminated square composite plates with SSSS (S2) 
boundary conditions subjected to in-plane uni-axial 
compression , plate thickness ratio (a/h)=20. (T*cri = 
Tcrie+08) - mean dimensionalized linear solution. It is 
seen that for given environmental conditions the 
increase in amplitude ratios increases the mean 
values and decreases the COV. The combined effects 
of all random input variables increase the COV, 
when plate is resting on Winkler nonlinear elastic 
foundations. On change of foundation parameters to 
Pasternak these values further increase and COV 
decreases. With the increase of moisture contents 
keeping temperature constant the effects of 
amplitude ratios decreases the mean values for 
Winkler nonlinear elastic foundations, however 
these values are higher for Pasternak nonlinear 
elastic foundations. Increase in environmental 
conditions drastically decreases the mean values and 
coefficient of variations. 



TABLE 6 (A) 







[0°/90°/]2T 


[0°/90°/]4T 






(Tcrri) 




COV, Xcml 




(Tcml) 




COV, Xcml 




Foundation 


Wmax/h 










































parameters 








bi 








bi 










(i=l..ll) 


(i-6,7) 


(i=8) 


(i=9..11) 




(i=l..ll) 


(i=6,7) 


(i=8) 


(1=9.-11) 


(Jfci=100, 


0.2 


0.0672 


0.2137 


0.0311 


0.1781 


0.0185 


0.0930 


0.1655 


0.0225 


0.1288 


0.0088 


ki=00 ,fo=100) 


0.4 


0.0859 


0.1781 


0.0243 


0.1393 


0.0089 


0.1218 


0.1465 


0.0171 


0.0983 


0.0136 




0.6 


0.1000 


0.1642 


0.0209 


0.1198 


0.0109 


0.1487 


0.1437 


0.0140 


0.0805 


0.0211 




Tcrl 


(0.0477) 










(0.0511) 










(fo=100, 


0.2 


0.0917 


0.1815 


0.0228 


0.1306 


0.0241 


0.1131 


0.1558 


0.0185 


0.1058 


0.0182 


fc=10, fe=100) 


0.4 


0.1060 


0.1642 


0.0197 


0.1129 


0.0191 


0.1422 


0.1438 


0.0147 


0.0842 


0.0178 




0.6 


0.1204 


0.1554 


0.0174 


0.0995 


0.0181 


0.1701 


0.1431 


0.0123 


0.0704 


0.0220 




Tcrl 


(0.0772) 










(0.0806) 










(fa=100, 


0.2 


0.0926 


0.1808 


0.0225 


0.1293 


0.0231 


0.1156 


0.1571 


0.0181 


0.1036 


0.0188 


fc=10, fo=200) 


0.4 


0.1652 


0.1000 


0.0193 


0.1105 


0.0198 


0.1498 


0.1501 


0.0139 


0.0799 


0.0236 




0.6 


0.1250 


0.1585 


0.0167 


0.0958 


0.0211 


0.1930 


0.1439 


0.0108 


0.0620 


0.0265 




Tcrl 


(0.0772) 










(0.0806) 
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TABLE 6 (B) 







[±45°]2T 


[±45°]4T 






(Tcml) 




COV, ^.crrj 




(Tcml) 




COV, Xcnl 




Hrtiitl/1ali/*in 
rOUIlLldliUIl 


V V max/ n 










































parameters 








bi 








bi 










(i=l..ll) 


(i-6,7) 


(i=8) 


(i=9..11) 




(i=l..ll) 


(i-6,7) 


(i=8) 


(i=9..11) 


(h=100, 


0.2 


0.0772 


0.1978 


U.Uz/U 


U.lCOU 


0.0066 


0.0996 


0.1617 


U.UzlU 




0.0038 


ki=00 ,fo=100) 


0.4 


0.0912 


0.1744 


0.0229 


0.1313 


0.0039 


0.1273 


0.1474 


0.0164 


0.0940 


0.0073 




0.6 


0.1051 


0.1626 


0.0199 


0.1140 


0.0054 


0.1499 


0.1473 


0.0139 


0.0799 


0.0121 






(0.0635) 










(0.0675) 










(h=100, 


0.2 


0.0894 


0.1907 


0.0234 


0.1339 


0.0128 


0.1106 


0.1633 


0.0189 


0.1082 


0.0105 


fa=io, fo=ioo) 


0.4 


0.1009 


0.1752 


0.0207 


0.1187 


0.0115 


0.1411 


0.1490 


0.0148 


0.0849 


0.0099 




0.6 


0.1160 


0.1620 


0.0180 


0.1032 


0.0104 


0.1662 


0.1491 


0.0126 


0.0720 


0.0127 




Tcrl 


(0.0809) 










(0.0849) 










(h=100, 


0.2 


0.0897 


0.1911 


0.0233 


0.1335 


0.0128 


0.1117 


0.1654 


0.0187 


0.1072 


0.0108 


fo=10, fa=200) 


0.4 


0.1017 


0.1770 


0.0205 


0.1177 


0.0118 


0.1450 


0.1564 


0.0144 


0.0826 


0.0132 




0.6 


0.1180 


0.1659 


0.0177 


0.1015 


0.0117 


0.1747 


0.1625 


0.0120 


0.0686 


0.0196 




Tcrl 


(0.0809) 










(0.0849) 











TABLE 7 (A) 









AT= 


300K AC=0.001 




AT= 300K AC=0.0025 






(Tcml) 




COV, Xcml 




(Tcml) 




COV, Xcml 




Foundation 


Wmax/h 










































parameters 








bi 








bi 










(i=l..ll) 


(i-6,7) 


(i=8) 


(i=9..11) 




(i=l..ll) 


(i-6,7) 


(i=8) 


(i=9..11) 


(fa=100, 


0.2 


2.3511 


0.2502 


0.0081 


0.0045 


0.0105 


1.0476 


0.2453 


0.0073 


0.0103 


0.0103 


fa=00 ,fo=100) 


0.4 


2.5527 


0.2512 


0.0074 


0.0041 


0.0093 


1.1319 


0.2458 


0.0068 


0.0095 


0.0121 




0.6 


2.7029 


0.2551 


0.0070 


0.0039 


0.0095 


1.1800 


0.2772 


0.0065 


0.0091 


0.0059 




Tcrl 


(2.1164) 










(0.8604) 










(fa=100, fa=10, 


0.2 


3.2970 


0.2367 


0.0058 


0.0032 


0.0176 


1.1909 


0.2372 


0.0064 


0.0090 


0.0200 


fe=100) 


0.4 


3.2297 


0.2369 


0.0059 


0.0030 


0.0171 


1.7236 


0.2249 


0.0044 


0.0062 


0.0098 




0.6 


3.1101 


0.2517 


0.0070 


0.0034 


0.0180 


1.2950 


0.2462 


0.0059 


0.0083 


0.0143 




Tcrl 


( 2.6525) 










(1.0782) 










(fa=100, fa=10, 


0.2 


2.8909 


0.2397 


0.0066 


0.0040 


0.0186 


1.1520 


0.2587 


0.0067 


0.0093 


0.0156 


fe=200) 


0.4 


3.8574 


0.2386 


0.0049 


0.0027 


0.0125 


1.6026 


0.2526 


0.0048 


0.0067 


0.0104 




0.6 


3.1339 


0.2385 


0.0061 


0.0036 


0.0176 


1.3021 


0.2468 


0.0059 


0.0083 


0.0143 




Tcrl 


(2.6525) 










(1.0782) 
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TABLE 7 (B) 







AT=325K AC=0.0125 


AT= 325K AC=0.015 










COV, A.cml 




^ 1 crnl f 




COV, Xcml 




Foundation 
parameters 


Wmax/h 






















(Tcrnl) 


bi 




bi 








Ci=1 1 1 1 


Ci=6 71 




fi=9 111 




Ci=1 1 1 1 


Ci=6 71 




fi=9 1 1 1 


(fa=100, 


0.2 


0.2197 


0.1342 


0.0076 


0.0545 


0.0037 


0.1832 


0.1390 


0.0076 


0.0654 


0.0037 


Kl— UU ,rij— 


0.4 


0.2543 


0.1286 


0.0066 


0.0471 


0.0040 


0.2120 


0.1324 


0.0066 


0.0565 


0.0040 




0.6 


0.2855 


0.1303 


0.0059 


0.0420 


0.0058 


0.2381 


0.1332 


0.0059 


0.0503 


0.0058 




Tcrl 


( 0.1906) 










(0.1589) 










(h=100, 


0.2 


0.2830 


0.1548 


0.0059 


0.0423 


0.0182 


0.2359 


0.1568 


0.0059 


0.0508 


0.0180 




0.4 


0.1490 


0.1000 


0.0065 


0.0463 


0.0155 


0.2842 


0.1362 


0.0049 


0.0422 


0.0116 




0.6 


0.2967 


0.1594 


0.0056 


0.0404 


0.0175 


0.2159 


0.1529 


0.0065 


0.0555 


0.0155 




Tcrl 


(0.2387) 










(0.1990) 










(fa=100, 


0.2 


0.2856 


0.1558 


0.0059 


0.0419 


0.0176 


0.2380 


0.1573 


0.0059 


0.0503 


0.0173 


fa=10, fa=200) 


0.4 


0.2704 


0.1451 


0.0062 


0.0443 


0.0146 


0.2882 


0.1357 


0.0048 


0.0416 


0.0113 




0.6 


0.3406 


0.1393 


0.0049 


0.0352 


0.0126 


0.4070 


0.1130 


0.0034 


0.0294 


0.0104 




Tcrl 


(0.2387) 










(0.1990) 











Conclusions 

The stochastic DISFEM procedure outlined in the 
present study has been used to obtain the mean and 
COV of the hygrothermal postbuckling load of the 
laminated composite plates without foundations and 
resting on nonlinear elastic foundations, subjected to 
uniform temperature and moisture rise with random 
system parameters. The following conclusion can be 
drawn from this limited study: 

1. For various input random variables, amplitude 
ratios of cross ply laminate, the increase in 
temperature and moisture concentration effect 
significantly as the hygrothermal postbuckling 
load decreases drastically specially for clamped 
support conditions. 

2. The random system properties, En, E22, G12, an, h 
and foundation parameters have dominant effect 
on the coefficient of variation of the hygrothermal 
post buckling load as compared to other system 
properties subjected to both temperature and 
moisture dependent hygrothermo-elastic 
properties. The strict control of these random 
parameters is therefore, required if high reliability 
of the laminated composite is desired. 



3. When the combined effects of aspect ratio and 
amplitude ratio along with random variables in 
hygrothermal environment are studied, then the 
clamped support conditions are significantly 
affected by the hygrothermal post buckling load 
compared to simple support conditions. However, 
coefficient of variation of hygrothermal 
postbuckling load for composite plates subjected 
to clamped support condition is lower as 
compared to simply support conditions for both 
without foundation and with nonlinear elastic 
foundations. 

4. The mean hygrothermal post buckling load for 
the parameters considered as mentioned above is 
more pronounced for square composite plates 
compared to rectangular plates; however it is vice 
versa in the case of coefficient of variation of 
hygrothermal postbuckling load for composite 
plates. The effect of hygrothermal load is more 
significant for coefficient of variation when 
number of layers increase. 

5. In general, the hygrothermal effects are of due 
importance for analyzing the structural response 
made of composite laminate. The negligence of 
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hygrothermal effects in analyzing the system 
behavior may leave the design unsafe for 
manufacturing and may prove to be a tool to 
avoid failure for reliable operational requirements. 
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